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Abstract 

A transition to Kraichnan ultimate regime of convection has been reported in very high Rayleigh 
numbers experiments, but not in all of them. These apparently contradictory results can be ex- 
plained by a recent phenomenological model which accounts for the non-ideality of the plate thermal 
properties [Chilld et al, Physics of Fluids 16:2452 (2004)]- In this paper, we present a direct test 
of this model, using a low conductivity plate. We found an unaltered transition, not compatible 
with the modeVs predictions. 



1 Motivation 

In 1962, Robert Kraichnan predicted 1 that a fully-turbulent regime of convection will settle in 
Rayleigh-Benard cells for an intense thermal driving force (very high Rayleigh number Ra). This 
regime is characterized by turbulent boundary layers and consequently an enhanced heat transport 
efficiency (Nusselt number Nu) compared to the other convection regimes. Kraichnan predicted a 

Nu ~ Ra^^'^ dependence, with a logarithmic correction which reduces the effective exponent right 
above the transition to this so-called ultimate regime. 

Only 5 experiments and one numerical simulation have explored Rayleigh numbers at least 
at Ra = 10^^ and higher. For all of them, the Prandtl number was in the intermediate region 
(0.7 < Pr < 20). The experiments have been conducted nearly in the same conditions : cylindrical 
cells of aspect ratio 0.5 or 1, cryogenic helium as the working fluid, copper top and bottom plates 
and thin stainless steel side wall. These very high Rayleigh number studies can be separated 
into two sets according to the measured Nu{Ra) dependence. For two experiments [SlEl, this 
dependence is smooth and follows approximatively a 0.3 exponent power law up to Ra ~ lO^'' 
for one of them 3 , while for two other experiments 01 [SI and the simulation 6 , the exponent 
increases from 0.3 to approximately 0.4 in the Ra = 10^^ — 10^^ region. This second regime has 
been interpreted as Kraichnan ultimate regime and this interpretation was strengthened by the 
measurement of the full 0.5 exponent in a rough surfaces cell j7|. 

No obvious parameter -such as the Prandtl number or the cell aspect ratio- seems to be 
correlated with the occurrence of the transition among these few studies. Therefore, understanding 
what can favor or suppress occurence of the ultimate regime remains today a major challenge for 
the understanding of very high Rayleigh number convection |H1 Ej • 

Recently, Chilla and colleagues proposed a phenomenological model ^Uj trying to explain 
these apparently contradictory results. This model asseses the non-ideality of the plates' thermal 
properties : if the conductivity and/or the heat capacity of the plates are not high enough, a 



constant temperature cannot be sustained at the fluid-plate interface during the development of 
thermal plumes. As a consequence, the global heat transfer is shown to saturate below Kraichnan 
prediction. The Chilla-rastello condition (Cr) to observe the transition is : 



Cr > 0.8 (1) 

Where Cr is a plate quality factor which depends on the ratios of plate to fluid conductivities 

heights ^ (a is the plate thickness) and heat capactities, as well as the usual dimensionless 

parameters Nu, Pr, aspect ratio F and Reynolds number Re. For the cryogenic experiments 
previously mentionned, a thin plate and quasistatic approximations can be applied to the general 
expression of Cr. One gets Jili : 

- h Xf RePr 

For reference, we should mention two other recent studies of flnite plate conductivity effect 
on convection ^iJEl^ although they are not directly related to the ultimate regime issue. Some 
previous works are reported in these two papers. 

In this paper, we present two experiments that have been conducted to test Chilla et al. model. 
We took out of the shelves two 20-cm high Rayleigh-Benard cells in which the transition have been 
observed in the past: one with smooth surfaces ^ and the other one with rough (corrugated) 
surfaces 7 . Both bottom plates were removed and replaced with two others, geometrically- 
similar but made out of brass: one with a smooth surface and the other with a rough one. The 
corresponding Cr parameter for these new bottom plates is smaller than 0.1 for Ra > 10^^. 
According to Chilla et al, the global heat transfer Nu{Ra) in the ultimate regime should be 
signiflcantly reduced compared to the original all-copper experiments. 



2 Experimental set-up and results 

In order to make the comparison with the all-copper experiments straightforward, the experimental 
set-up and the operating procedure are very similar to those of these previous studies HI • 
Since this experimental information is detailed in the previous papers, we shall only report here 
the differences with these experiments. 

The conduction of the brass plates has been measured with the inductive technique already 
used to measure the conductivity of the copper plates [l3j. We found ksrass = 4:.0Wm~^K~^ 
±10% at 4.2 K, to be compared with kcu = 1090 Wm-^i^^"^ and 880Wm-^K-^ at 4.2 K for 
the smooth surfaces and rough surfaces copper plates. In our experimental conditions, the heat 
capacity of brass is comparable to the one of copper ^Sl • The surface roughness of the corrugated 
plate has been measured with a profllometer. The grooves are 450 /im apart JHl? 145 /im deep 
and look like the reference copper plate roughness (see flg. Both 25mm-thick brass bottom 
plates are heated from below. In order to minimize the heater footprint on the other side of the 
plates, the heating wire has been uniformly and carefully varnished with a 1.5 mm spacing on the 
bottom surface. 

A special attention has been dedicated to the temperature measurements on the axis of the 
brass plates, near the fluid-brass interface. Finite elements simulations ^\ have been used to 
optimize the shape of the copper thermometer holders. Optimization was performed for 10^ < 
Ra < 10^^ and by modeling the fluid's thermal boundary layer as a diffusive and spatially uniform 
slab of thickness h/2Nu, where h is the height of the cell. The thermometer holders are inserted 
into the brass plates, in long 04 mm holes ending with a 3 mm long 03 mm thread, itself ending 
2 mm away from the surface in contact with the fluid. The holes and the thermometer holders 
are responsible for a local non-uniformity of temperature at the fluid-brass interface. According 
to the simulation, and for all the data presented in this paper, the magnitude of this effect is 
less than 1% of the temperature difference across the fluid. To this extend, the temperature 
measurement can be assumed as non-invasive. A correction needs to be applied to compensate for 
the temperature drop between the thermometer holders and the fluid-brass interface. According to 
the simulation, this correction represents a 0.11 K/W resistance, with negligible Nu dependence, 
which corresponds to a effective position of the probe 3.5mm from the fluid-brass interface. For 
all the data presented in this paper, this correction represents less than 20% of the temperature 
difference across the fluid (less than 10% for the smooth surfaces cell). 

Temperature on both the top and bottom plates have been measured with doped Ge resistors, 
in-situ calibrated : 1) with respect to each other, 2) with respect to a commercially calibrated 
one |18j, 3) with respect to the fluid critical point. The estimated uncertainty on the temperature 
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Figure 1: Profilometer measurement of the roughness profile along a bottom plate diameter. 
Distances are given in microns. Upper plot : copper plate, Lower plot : brass plate. 

difference between the top and bottom plate is 1 mK. All the data presented have been obtained 
with a temperature difference across the fluid larger than 20 mK. 

The arbitrarily Boussinesq criterion is that no fluid parameter should var y by more than 20% 
across cell. Recent improvements in the knowledge of helium fluid properties jl9j have been taken 
into account in the calculation of Ra and Nu. Nu has been corrected for the side wall conduction 
using the analytical correction formula presented in 20 . The magnitude of this correction is 
nearly the same for all the cells. It decreases very rapidly with Ra and represents less than a 5% 
correction above Ra = 10^^. 

Figure [21 and [2l shows the heat transfer data in the rough and smooth surfaces "brass" experi- 
ments. The various solid symbols are associated with data acquired for a fixed cell density and a 
nearly fixed average temperature in the cell. The Prandtl number associated with each data set 
is in the 0.7 < Pr < 5.4 range. For comparison, the all-copper experiment data, restricted to the 
same Pr window, are also plotted (open circles). 

On Figure [3 the rough surfaces Nu data are compensated by Ra^^'^ so that Kraichnan's 
Nu{Ra) power law corresponds to a plateau. Such plateaux are found for both the brass and the 
all-copper experiments, and -within experimental uncertainty- they fall on each other. Some extra 
data between Ra = 3.10^^ and Ra = 6.10"^^ are not shown because they don't meet our various 
criteria (Boussinesq, Prandtl window, temperature-drop maximum correction and minimum tem- 
perature difference across the fluid). Nevertheless these extra data for the two experiments keep 
falling on each other. 

On Figure 13 the smooth surfaces Nu data are now compensated by Ra^^^. To our knowledge, 
the 1/3 exponent is the highest exponent predicted by the theories of turbulent convection, let 
aside Kraichnan's regime. Above Ra = 4.10^^, if we fit our data with a local power law, we find 
an effective exponent of dlog{Nu) / dlog{Ra) ~ 0.4, significantly larger than the 1/3 prediction. 
The scattering of the data from the reference all-copper experiment is larger because the temper- 
ature difference across the cell hasn't been average as long. Here again, within the experimental 
uncertainty, the data from both the brass and the all-copper experiments fall on each other. 



3 Concluding remarks 

The two experiments presented here are intended to test Chilla et al. model ^j, which was able 
to explain why the ultimate regime is observed in some experiments but not all of them. Within 
experimental uncertainties, we find that the substitution of a copper plate by its replica made out 
of brass has no consequence on the global heat transfer, in particular within the ultimate regime. 
This result is in clear contradiction with the model. 

It seems difficult to adjust the model parameters while preserving its ability to discriminate 
between the cryogenic experiments |2l 0] |2l Hi Ej • This suggests that the critical difference between 
these experiments has to be searched elsewhere. A candidate is the large scale flow structure 
(LSFS). Its bimodality, which was evidenced [2IE21 foi" < 10^^ i^i^y extend to higher Ra as 
suggested by the two figures of this paper. In 2001, it was conjectured 14 that this transition to 
the ultimate regime could only occur for one LSFS, at least below Ra ~ 10^''. 




Figure 2: Compensated Nusselt number Nu/Ra^^'^ versus Ra for the rough surfaces experiments. 
Open circles : copper bottom plate, Solid symbols : brass bottom plate. 




Figure 3: Compensated Nusselt number Nu/Ra^^^ versus Ra for the smooth surfaces experiments. 
Open circles : copper bottom plate, Solid symbols : brass bottom plate. 
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